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Изложена численная методика рационального прочностного проектирования тел вращения при осесимметричном нагружении. 
Результатом такого проектирования являются конструкции, которые в литературе принято называть дискретно равнопрочными. 
Целевая функция как таковая при этом не минимизируется, ее значения лишь иллюстрируют ход итерационного процесса, т.е. 
оптимизационная задача заменяется параметрическими исследованиями. Приводятся примеры прочностного проектирования 
метаемых элементов для высокоскоростных баллистических установок в квазистатической постановке. 


1. Введение 

Проблемы оптимизации конструкции, имею- 
щие большое прикладное и теоретическое значе- 
ние, привлекают большое внимание исследовате- 
лей. Из работ, посвященных вопросам оптимиза- 
ции массы и формы упругих тел и тонкостенных 
конструкций с позиций аналитического и числен- 
ного решения задач математического программи- 
рования, следует отметить монографии [1-3]. 
Уменьшение массы конструкции без потери ее 
прочностных качеств является своего рода "сверх- 
задачей" для любого проектировщика, если не ос- 
новной целью. Эта цель достигается путем выбора 
рациональной схемы конструкции с необходи- 
мостью, кроме других ограничений, выполнения 
требований по прочности. Решение задач опти- 
мального проектирования в механике деформируе- 
мого твердого тела (МДТТ) приводит к сложным 
математическим задачам и наталкивается на серь- 
езные математические трудности. Смысл построе- 
ния и использования упрощенных методик сводит- 
ся к тому, что они значительно сужают исследуе- 
мую область изменения управляющих параметров, 
в которой находится оптимальный проект, а также 
позволяют получать начальное приближение для 
использования методов математического програм- 
мирования или экспериментальной отработки 
проектов конструкций. 

В литературе [2] термины "оптимальное" и "ра- 
циональное" проектирование принято различать 
следующим образом: если задача решается в стро- 
гой математической постановке, то говорят об оп- 
тимальном проектировании, если же речь идет о 
некоторой приближенной процедуре, это принято 
называть рациональным проектированием. Под 
второе определение подходит описанная ниже, а 
также в [4, 5], процедура отыскания лучшего в 
смысле прочности проекта тела вращения. 

2. Алгоритм рационального прочностного 

проектирования тела вращения 

Задача о минимуме массы конструкции при вы- 
полнении условий прочности в той или иной фор- 
ме, также как и рациональное распределение мате- 


риала заданного объема, имеют большое практи- 
ческое значение при создании различных 
конструкций. 

В качестве параметра, характеризующего нап- 
ряженно-деформированное состояние (НДС) 
конструкции и эффективность работы материала в 
отдельных точках, принимается интенсивность де- 
формаций е,. Задача формулируется следующим 
образом: найти такое распределение материала по 
конструкции, при котором конструкция имеет ми- 
нимальную массу, а величина е, во всех точках не 
превосходит заданное значение е тах . Для осесим- 
метричных конструкций, которые можно рассмат- 
ривать как толстостенные оболочки вращения, за 
управляющий параметр принимается толщина И 

к(г) = /?і (2)-К 2 (2), 

которая при фиксированной отсчетной - внешней, 
внутренней или срединной - поверхности задает 
распределение материала по конструкции. и 
К 2 (і) представляют собой радиальные координаты 
наружной и внутренней поверхности соответствен- 
но. В дальнейшем будем называть функцию Н(і) 
толщиной. При осесимметричном нагружении те- 
ла вращения внешними нагрузками задача сводит- 
ся к определению функции И(і), дающей минимум 
функционалу 




Ѵ = тг^[К^(г)-К;(г)]сІ2 

0 

(1) 

при ограничениях прочностного и геометрическо- 

го характера 


е < е 

і тах ’ 

(2) 

/ 7 ( 2 )>/ 7 0 . 

(3) 


Здесь V - объем тела вращения; I - размер в 
осевом направлении; е тах - максимальное допусти- 
мое значение интенсивности е;, И 0 - минимальное 
допустимое значение толщины, назначаемое про- 
ектировщиком. 

Соотношения (1-3) вместе с уравнениями рав- 
новесия, условиями на поверхности, геометричес- 
кими и физическими соотношениями представля- 
ют собой формулировку задачи нелинейного мате- 
матического программирования. 
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В сформулированной задаче проектирования 
нет явно выраженной связи между значениями уп- 
равляющей функции и ограничениями прочност- 
ного характера (2). В связи с этим задача перефор- 
мулируется следующим образом. При численном 
решении задачи определения НДС на осевое сече- 
ние конструкции наносится двумерная конечно- 
разностная сетка, разбивающая тело по длине на 
так называемые подконструкции. Тем самым уп- 
равляющая функция к(і), которую можно тракто- 
вать как вектор бесконечной размерности, заменя- 
ется набором величин /г,=/г(г,), каждая из которых 
постоянна в пределах /-ой подконструкции (/ = 1, 
2 ,...,/; / -число разбиений тела по длине) и имеет 
смысл толщины. Этой заменой задача приводится 
к управлению вектором конечной размерности /. 
Распределение толщин /г, находится методом пос- 
ледовательных приближений по формуле 

Ь- к+ 4 = тах[/? 0 , /гД \\-Ъ + Ъе^ / е тах )], 

(4) 

г-1,2, А = 0,1,2,..., 

где к - номер итерации; е іти - максимальная ин- 
тенсивность деформаций в /-ой подконструкции, 
разбиваемой в свою очередь на К частей по толщи- 
не; Ь - своего рода коэффициент релаксации, вве- 
дение которого позволяет ускорить сходимость 
итерационного процесса. 

Подобный способ определения переменной оп- 
тимальной толщины для тонкостенных конструкций 
по оболочечной теории использован в работах [6-8]. 

За нулевое приближение берется некоторая 
постоянная толщина Н. Если вновь найденная тол- 
щина /г) Ш) не удовлетворяет условию (3), то она при- 
равнивается /г 0 . Формула (4) имеет следующий 
смысл: на каждом шаге расчета толщина /г, /-ой под- 
конструкции увеличивается, если интенсивность 
деформаций больше заданной е тах , в противном слу- 
чае уменьшается, если не нарушается условие (3). 

Таким образом, алгоритм прочностного проекти- 
рования представляет собой многошаговый итераци- 
онный процесс. На каждом шаге определяются пара- 
метры НДС, в том числе интенсивность деформаций 
е„ и затем проводится перераспределение материала 
по конструкции по формуле (4) с учетом ограниче- 
ний (2), (3). Процесс вычислений заканчивается, ес- 
ли одновременно выполняются соотношения 

|Г<* + 4 - ѵ™ ) / Г к+ 4 1 < и | (еМ і - с тах )/е тш \<е е (5) 

или условие 

\(Н ср ~К)! ^о| < е іі > (6) 

где Ѵ (к) , Ѵ (к+1> - объемы тела вращения, соответству- 
ющие распределениям толщин к (к> , /г (Ш) ; е ѵ , е е , г,, - 
погрешности определения объема и выполнения 
ограничений (2) и (3) соответственно; еМ, - макси- 
мальная величина из I величин е /тах ; Н ср - средняя 
толщина тела вращения, определяемая как среднее 
арифметическое толщин I подконструкций. Пос- 
леднее условие необходимо для окончания итера- 
ционного процесса при выходе значений толщин 


Ні на ограничение А 0 , что говорит об излишнем за- 
пасе по толщине, задаваемом условием (3). 

Результатом такого проектирования являются 
конструкции, которые в литературе принято назы- 
вать дискретно равнопрочными. Целевая функция 
как таковая при этом не минимизируется, ее значе- 
ния лишь иллюстрируют ход итерационного про- 
цесса, т.е. оптимизационная задача заменяется па- 
раметрическими исследованиями. 

Определение осесимметричного упругопласти- 
ческого НДС проводится вариационно-разност- 
ным методом (ВРМ) [9, 10] с использованием 
быстродействующего метода Гаусса [11-13] для ре- 
ализации системы линейных алгебраических урав- 
нений большого порядка относительно искомых 
перемещений. Физически нелинейная задача ре- 
шается с помощью соотношений деформационной 
теории пластичности и метода переменных пара- 
метров упругости. 

3. Примеры проектирования конструкций 

при высокоскоростном метании 

Изложенный алгоритм прочностного проекти- 
рования иллюстрируется примерами квазистати- 
ческого расчета осесимметричных конструкций 
при высокоскоростном метании. 

С точки зрения МДТТ процесс разгона 
конструкции является достаточно длинным во вре- 
мени - несколько миллисекунд, и из расчетной 
практики известна возможность квазистатического 
подхода к оценке прочности большинства подоб- 
ных конструкций, что связано с согласованием ха- 
рактерного времени нагружения и периода пробега 
упругих возмущений по конструкции малых габа- 
ритов. Решение исходной динамической задачи в 
квазистатической постановке возможно благодаря 
введению в исходную расчетную схему нагружения 
массовых сил, являющихся следствием использо- 
вания связанной с метаемым телом неинерциаль- 
ной системы координат [14-16]. Для работы с по- 
являющимися фиктивными массовыми силами 
вводится понятие коэффициента перегрузки (или 
просто перегрузки /V), который определяется отно- 
шением ускорения тела а в данный момент к уско- 
рению силы тяжести §: /Ѵ= а/§. 



Рис. 1 . Конфигурация исходного и оптимального (заштри- 
хованная область) проектов 

Рассматривается (рис. 1) тело вращения (мате- 
риал - Ст. 3) постоянной толщины Н=АВ= 0,5 см, 
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срединная поверхность которого образована 
частью эллипсоида вращения с полуосями 5,0 см и 
37,6 см. Длина I =30,0 см. Тело нагружено массо- 
выми силами с ІѴ=— 5000. За отсчетную принимает- 
ся срединная поверхность тела. Вследствие сим- 
метрии конструкций относительно оси вращения 
0г на рисунках представлены лишь половины осе- 
вых сечений рассматриваемых тел. 



Рис. 2. Распределение максимальных интенсивностей де- 
формаций е^ах в 1) исходном и 2) оптимальном про- 
ектах подлине тела 

Конфигурация исходного и оптимального (зашт- 
рихованная область) в указанном выше смысле про- 
ектов для ограничений А 0 = 0,1 см, е ти =0,002 (что со- 
ответствует деформации начала текучести) представ- 
лены на рис. 1. Характер изменения максимальных 
интенсивностей деформаций е іт для обоих проек- 
тов показан на рис. 2. В отличие от исходного (кри- 
вая 1), в полученном оптимальном проекте (кривая 
2), согласно ограничению (2), материал деформиру- 
ется упруго. Радиальные перемещения внешней по- 
верхности положительны и их наибольшие значения 
у торца г=0 достигают 0,00834 см для исходного и 
0,004 см для оптимального проекта. Максимальные 
осевые перемещения торца і=Ь отрицательны и 
имеют соответственно значения 0,077 см и 0,053 см. 
Объем оболочки исходного проекта конструкции с 
417 см 3 сократился на 67 % до 138,09 см 3 , тем самым 
значительно увеличив полезный объем внутри тела 
вращения. Представленные результаты получены за 
десять итераций (под итерацией понимается расчет 
НДС тела с распределением толщины, полученным 
по формуле (4) для гк=г г =е,,=0,01. 



Рис. 3. Конфигурация исходного и оптимального (заштри- 
хованная область) проектов 

На рис. 3, 4 представлены результаты проектиро- 
вания толстостенного тела вращения для случая при- 
ложения массовых сил и давления. Геометрия внеш- 
него контура определялась как часть усеченного эл- 
липса с полуосями 8,0 см и 38,35 см; длина равна 
30,0 см, толщина исходного варианта конструкции 


#=2,0 см, перегрузка #=-3000, амплитуда приложен- 
ного на части внешней поверхности давления 
Р=-500 МПа. Материал корпуса - сталь с деформаци- 
ей начала текучести ^=0,003; допустимая деформация 
^ітж = 0,05; #=0,5 см. Здесь и в следующем примере за 
отсчетную принимается внешняя поверхность тела. 
На торце г=0, как и в предыдущем рассмотренном 
примере, выполняются условия жесткой заделки. Ко- 
нечное распределение толщины (заштрихованная об- 
ласть на рис. 3) получено за пять итераций. 



Рис. 4. Распределение максимальных интенсивностей де- 
формаций е ішш в 1) исходном и 2) оптимальном про- 
ектах по длине тела 

Масса тела в конечном варианте составляет 32 % 
от массы исходного проекта. Такое уменьшение 
стало результатом "срезания" материала в менее 
нагруженных частях конструкции и значительного 
"наращивания" материала (аналог силового элемен- 
та типа шпангоута) в месте приложения внешнего 
давления. Имеет место значительное увеличение 
интенсивности деформаций материала в опти- 
мальном проекте, величина которой в районе этого 
"шпангоута" вышла с заданной погрешностью 
е с =0,05 на ограничение е тх . Пунктирная линия на 
рис. 4 обозначает деформацию начала текучести е 5 . 

В качестве последнего примера применения 
представленного метода рационального прочност- 
ного проектирования рассмотрена модельная зада- 
ча о метании поддона сложной формы. Для иссле- 
дования степени эрозии новых конструкционных 
материалов при движении в атмосфере с высокой 
скоростью используются установки, позволяющие 
с помощью специальных поддонов метать образцы 
материалов со скоростью 3...6 км/с. На рис. 5 
представлена расчетная схема выполненного из 
дюралюминия такого поддона. Исследуемый обра- 
зец материала помещается на внешней поверхнос- 
ти ОД дна поддона. Масса поддона равна 120 г, мас- 
са метаемого образца - 50 г. Исходная конфигура- 
ция подлежащего корректировке внутреннего кон- 
тура поддона обозначена ИНОЕ. Вычисления про- 
водились для следующей геометрии поддона: 
#=1,75 см; #=2,25 см; #=2,5 см; #=2,0 см; 
1=5,0 см; 0Р=4,0 см; 0М=0,7 см; ВС=ВЕ= 0,5 см. Ге- 
ометрические ограничения: для дна 0,5 см; для ци- 
линдрической части корпуса #=2,0 см; #=2,2 см. 

На внутреннем МНСР и внешнем ВСОЕ конту- 
рах поддона действует давление газа Р, посред- 
ством которого метается поддон. Величина перег- 
рузки, приложенной в противоположном движе- 
нию поддона направлении, вычислялась из усло- 
вия равенства нулю вектора всех действующих 
вдоль оси 0", сил, и для Р =100 МПа равнялась 
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7Ѵ= 11600. Давление на поверхность ОД от присое- 
диненной массы Р п ~ 39 МПа. На контактных с ме- 
тательной установкой поверхностях АВ и ЕР поддо- 
на ставятся условия отсутствие трения и ограниче- 
ния на радиальные перемещения. Особенностью 
постановки граничных условий для данной задачи 
является то, что на осевые перемещения не накла- 
дывается ограничений. 



Рис. 5. Расчетная схема метаемого поддона и результаты 
проектирования. Конфигурация внутреннего контура 
поддона: МНОР — исходная; МіфР — оптимальная; 
КІОР — оптимальная при отсутствии образца иссле- 
дуемого материала 

Расчет НДС исходного проекта выявил упругий 
характер работы материала метаемого поддона, осо- 
бенно в его цилиндрической части. Часть дна у оси 
вращения близка к выходу в пластическую область 
деформирования. Использование разработанной 
методики позволило получить новую конфигурацию 
МІО Г внутреннего контура поддона. На рис. 5 вид- 
но, что его размеры вышли на заданные геометри- 
ческие ограничения по всей его поверхности, за иск- 
лючением небольшой части дна на оси вращения 0?. 
Почти весь материал дна поддона деформируется 
пластически, но при этом интенсивности деформа- 
ций не превосходят заданное ограничение е тах =0,02. 
Вся цилиндрическая часть корпуса находится в об- 
ласти упругих деформаций, что указывает на запас 
по прочности для материала в этой части поддона. 

Анализ НДС позволяет сделать вывод о наиболее 
напряженной работе материала дна поддона. При- 
мерно две трети его внешней поверхности получают 
осевые перемещения по направлению метания под- 
дона, которые почти по линейному закону уменьша- 
ются от наибольшего значения ѵѵ=— 0,05 см на оси 
вращения О^до положительного значения тѵ=0,012 см 
в точке А. Величины радиальных перемещений цили- 
ндрической части корпуса поддона на внутренней 
поверхности имеют почти нулевые значения, а на по- 
верхности СП они отрицательны (т.е. направлены к 
оси вращения) и также малы. Интенсивности напря- 
жений примерно одинаковы как вдоль дна, так и по 
его толщине. На его внутренней поверхности их 
средняя величина равна —240 МПа. 

Итогом применения алгоритма рационального 
прочностного проектирования явилось уменьше- 
ние массы метаемой сборки на 24 %, а массы под- 
дона -на 34 %, которая в окончательном варианте 
составила 79 г. Величина перегрузки увеличилась 
на 32 % и достигла 7Ѵ=15200. 


На рис. 5 линией АУ (ГГ представлен внутренний 
контур оптимального проекта поддона, получен- 
ного для случая метания его без присоединенной 
массы. Такое изменение конфигурации внутрен- 
ней поверхности дна объясняется отсутствием при- 
ложенного на поверхности 0А давления и, как 
следствие этого, увеличением в материале дна зна- 
чений интенсивностей деформаций, что, согласно 
алгоритму проектирования, автоматически приво- 
дит к увеличению толщин подконструкций. Вели- 
чина перегрузки для этого варианта увеличивается 
с 16350 до 20200, а масса уменьшается до 84 г. 

4. Выводы 

Ввиду того, что при определении НДС система 
линейных алгебраических уравнений, являющаяся 
следствием использовании условия стационарнос- 
ти функционала полной потенциальной энергии 
системы, реализуется быстродействующим мето- 
дом Гаусса, представленный алгоритм рациональ- 
ного прочностного проектирования позволяет с 
относительно небольшими временными затратами 
для конкретного нагружения в статической и ква- 
зистатической постановках получать рекоменда- 
ции по конструктивному выполнению тел враще- 
ния при заданных ограничениях прочностного и 
геометрического характера. В рассмотренных при- 
мерах результаты были получены на относительно 
грубых сетках за несколько минут. 

Обладая достаточной наглядностью и простотой 
реализации, алгоритм позволяет получать проекты 
меньшей массы путем изменения толщины. Прове- 
денные расчеты для разных конструкций метаемых 
элементов и разных комбинаций их нагружения по- 
казали, что независимо от того, является НДС ис- 
ходной конструкции упругим или упругопластичес- 
ким, разработанный алгоритм приводит к уменьше- 
нию объема по сравнению с первоначальным вари- 
антом, если он взят с запасом по прочности. В спро- 
ектированной конструкции исключается недопус- 
тимое НДС при рассматриваемом уровне нагрузки, 
что является следствием выполнения ограничения 
(2), которое эквивалентно принципу равнопрочнос- 
ти. В точках, где ограничение (2) выполняется в ви- 
де неравенства, возникают так называемые "мерт- 
вые" зоны, "... что приводит ... к нерациональному 
использованию материала несущих конструкций" 
[2]. По-видимому, при проектировании невозможно 
создавать сложные конструкции без таких зон, не 
вступая в противоречие с требованиями непрочно- 
стного характера - технологичности, эффективнос- 
ти, аэродинамики, баллистики и т.д. 

Описанный алгоритм не является реализацией 
какого-либо из методов математического програм- 
мирования, приводящего к получению оптималь- 
ной в известном смысле конструкции. Тем не ме- 
нее, результаты, полученные с его помощью, дос- 
таточно наглядны, физически оправданны и могут 
помочь конструктору при создании образцов но- 
вой техники. 
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